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Abstract. In order to increase the corrected field of view of an adaptive optics (AO) system, 
several deformable mirrors (DM) have to be placed in the conjugate planes of the dominant 
turbulent layers (multi-conjugate adaptive optics, MCAO (Beckers, 1988)). The performance 
of MCAO systems depends on the quality of the wavefront sensing of the individual layers and 
on the number of corrected modes in each individual layer as in single layer AO systems. In 
addition, the increase in corrected field of view depends on the number of guide stars providing 
information about the turbulence over a sufficiently large area in each turbulent layer. In this 
paper, we investigate these points and provide formulae for calculating the increased field of 
view with a new approach using the spatial correlation functions of the applied polynomials 
(e.g. Zernike). We also present a new scheme of measuring the individual wavefront distor- 
tion of each of the dominant layers with a Shack-Hartmann-Curvature Sensor using gradient 
information as well as scintillation. An example for the performance of a two layer MCAO 
system is given for the 3.5-m telescope of the Calar Alto Observatory, Spain, using a measured 
C^-profile. The corrected field of view in K-band (2.2pm) can be as large as 3 arcmin with a 
Strehl ratio above 60%. 
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1. Introduction 

While AO systems increase the angular resolution of ground based telescopes 
by a factor of 5-20, their most severe disadvantage is the very small corrected 
field of view which is typically of the order of 30" in K-band (2.2 /am) and 
only a few arc seconds in the visible. MCAO systems with several DMs in 
the conjugate planes of the dominant turbulent layers as proposed by Beckers 
(1988) have two important advantages over conventional AO systems: 

— The corrected field of view is increased considerably by correcting the 
dominant turbulent layers instead of correcting the integrated wavefront 
aberrations even if the turbulence is not located in distinct layers. 
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— Focal anisoplanatism which is a severe limitation for Laser Guide Star 
(LGS) AO systems can be almost completely eliminated by using mul- 
tiple LGSs (see Sect. 2.2). 

The optimal position of the DMs is determined by the -profile and by 
the system characteristics e.g. the size of the corrected field of view and 
the number of corrected modes. Usually the angular anisoplanatism of AO 
systems is calculated according to the (6/60 )' 5 / 3 ) -law (Tyson, 1991). For 
real systems however the corrected field of view is larger than given by this 
formula. Therefore a more realistic calculation is needed, both for estimating 
the performance of AO and MCAO systems, i.e. the corrected FOV, and for 
obtaining the optimal position of the DMs in MCAO systems. This is done in 
section ||, where we derive an equation for calculating the anisoplanatism both 
for AO and MCAO systems using correlation functions of the applied Zernike 
polynomials. The anisoplanatism depends on the C} r -profile, the number of 
corrected modes of each layer, the number and position of the DMs and the 
angle between science object and central guide star. Using various MCAO 
geometries, we apply these formulae to calculate the Strehl ratios for a 3.5-m 
telescope and for a turbulence profile measured by Kliickers (Kliickers et al., 
1998) at the Calar Alto Observatory. 

The second important issue of MCAO is the problem of wavefront sensing 



distinguishing the turbulent layers. Section 3.1 gives a short introduction to 
scintillation and describes the principle of separating the wavefront distor- 
tions of two turbulent layers with a Shack-Hartmann (SH) sensor (Ribak, 
1996; Schwartz et al., 1994) or a Shack-Hartmann-Curvature (SHC) sen- 
sor (Glindemann and Berkefeld, 1996), both using intensity fluctuations. By 
varying the effective altitude of the wavefront sensor(s) the SNR of the scin- 



tillation signal can be optimized. It will be shown in Section that a SHC 
sensor allows for much fainter guide stars than a SH sensor. Furthermore, the 
SH sensor delivers only one scintillation signal, therefore wavefront modes 
with zero-Laplacian cannot be properly measured. With a SHC sensor, how- 
ever, the second scintillation signal can reduce the error induced by these 
modes. 

A good introduction for the concept of MCAO can be found in Beckers 
(1988), more recent investigations of its applicability and performance were 
presented by Johnston and Welsh (1994). 

Unless mentioned otherwise, the wavelength used corresponds to the K- 
band (2.2/um). 
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2. The wavefront error a\ of angular anisoplanatism 

The wavefront error due to angular anisoplanatism is denoted by a 2 ,. It is 
caused by slightly different turbulence patterns in the light paths of two stars. 
Gg can be expressed as 



where 6 denotes the angular separation of the two stars. 6o is the isoplanatic 
angle that is defined such that a 2 , < 1 for objects lying inside the isoplanatic 
patch (Tyson, 1991). 6o has to be defined in such a way that the number of cor- 
rected modes is taken into account: A simple tip/tilt system has a much larger 
isoplanatic angle than higher order AO-systems. For the special cases of a 
tip/tilt system and certain higher order systems Sasiela and Shelton (1993) 
and Chassat (1989) already derived expressions. 

2.1. Angular Anisoplanatism of AO-systems 

The following approach of calculating o 2 , is based on the cross-correlation 
function C Xj y of two tip/tilt measurements at b and at b' and includes the 
general dependence of the isoplanatic patch on the degree of correction. 

Measurements of correlation functions yield important information about 
the parameters of AO systems: The temporal autocorrelation of tilt mea- 
surements determines the necessary bandwidth of the control loop and the 
correlation of tilt measurements of different stars (spatial correlation) allows 
us to calculate of the isoplanatic angle. 

The differential jitter of two beams caused by anisoplanatism is (Valley 
and Wandzura, 1979) 

((b-b')% y = 2[l-C X;y ]-(b 2 ), (2) 

where x- and y- directions correspond to the directions parallel and perpen- 
dicular to the line of sight between the two stars, (b 2 ) represents the one 
dimensional variance of the tilt. The variance of the differential image motion 
can be written as 



(AG) 2 = {{b-b'f) x + {{b-b'f) y 

= ((l-C x (e)) + (l-C,(6)))-2<£ 2 >, (3) 

where 2(b 2 } is the two-axis variance of the position angle b. The variance 
(b 2 ) of the position angle is related linearly to the wavefront error due to G- 
tiltn. In a very similar fashion, one can write the wavefront error a 2 , tt of the 

1 The wavefront error related to the motion of the image centroid is often called G-tilt since 
it is equivalent to the average gradient of the wavefront over the telescope aperture. 
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Zernike tilt - the first mode (apart from piston) of the Zernike decomposition 
- due to differential image jitter as 



<tt = ((1-Ct,,x(e)) + (1-Ctt, y (0)))-2.O.448 ^-J 
/D\ 5/3 

= 2- - • 0.448 • (2 - C tt , x (6) - C tt , y (6) ) . (4) 

0. 448(D/ro) 5 / 3 is the wavefront variance of the single axis tilt as given by 
Noll (1976). As soon as the correlation functions C tt . x ,y (0) have a value smaller 
than 0.5 the resulting tip-tilt correction with an on axis guide star deteriorates 
the image quality since the wavefront error becomes larger than 0.448(D/Vo) 5//3 , 

1. e. larger than the uncorrected value. Eventually, for zero correlation the 
variance is twice as large as without correction. 

The anisoplanatism of higher order wavefront modes can be deduced in a 
similar fashion (see Berkefeld (1998), corrected for minor inaccuracies in the 
deduction). Then , the wavefront error Gg N as a function of the correlation of 
N corrected Zernike modes can be written as 

al N = 2 (£) 7 • £(<4j-i - <4j)(l -CM), (5) 

summing up the different correlations Cj(Q) of each mode j and the modal 
variance a| t ._j — g\ x • to the wavefront error. The total wavefront error for a 
star at a distance 6 from the guide star can then be calculated by adding Gq n 
to the fitting error Gy it j >N . Eq. || can be used with any set of polynomials 
if the correlation functions and the modal variances are available. Chassat 
(1989) presented results for Gq N for selected Zernike-polynomials. 

Table | shows a| t for Kolmogorov turbulence described by Zernike poly- 
nomials as a function of the number j of corrected modes (see Noll (1976)). 

For j > 20, (Jg t ■ can be approximated by 

Gl tJ «0,2944(j + l)-^/ 2 (D/ro) 5 / 3 . (6) 

For an AO system with an LGS at infinity, i.e. without the cone effect, and an 
NGS for tip-tilt the wavefront error Gq n lgs becomes 

/ n\ V 3 

oJUlgs = 2 • {-) [0.448 • (2 - C X (Q NGS ) - C y (Q NGS )) + 

-<;)(! -Cj(Qws))}. (7) 

;=3 
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Table I. Wavefront error for Kolmogorov-turbulence described by Zernike polynomi- 
als as a function of the number j of corrected modes (Noll (1976)). 





= 1.0299(Z)/r ) 5 / 3 


<7 


= 0.0525(£>/V ) 5/3 


2 

"fit, 14 


= 0.0279 (Z?/r ) 5/3 


<1 


= 0.5820(D/r ) 5 / 3 
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= 0.0463(£>/V ) 5/3 


C fit,15 


= 0.0267(£>/V ) 5 / 3 


<2 


= 0.1340(£>/r ) 5/3 
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= 0.0401 (£>/V()) 5/3 


°fit,16 


= 0.0255(D/r ) 5/3 


<3 


= 0.1110(£>/r ) 5/3 


a 2 
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= 0.0377(D/Vo) 5/3 


°ut,17 


= 0.0243(£>/r ) 5/3 
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= 0.0880(£»/r ) 5 / 3 


<n 


= 0.0352(£»/r ) 5 / 3 


a 2 

"fit, 18 


= 0.0232(£»/r ) 5 / 3 


<4,5 


= 0.0648(£>/r ) 5/3 


a 2 

a fit,12 


= 0.0328(D/r ) 5/3 


a 2 

"fit. 19 


= 0.0220(£>/r ) 5/3 


<6 


= 0.0587(£»/r ) 5 / 3 


C fit,13 


= 0.0304(D/r ) 5/3 







The correlation functions Cj(B) of the individual Zernike modes j are 



cm 



J C j(zQ)C 2 n (z)dz 
fCldz ' 



(8) 



where Cj(zB) = cj(d) denotes the correlation caused by a single layer at al- 
titude z (Valley and Wandzura (1979)). Figure [T] shows the correlation for 
different aberrations as a function of the parameter d/D. Higher order aber- 
rations generally decorrelate faster because they correspond to smaller struc- 
tures (and thus smaller correlation lengths). 



correlation 



ci (y-tilt) 




C7 (coma) 



Figure 1. Correlation functions Cj{d) for different Zernike-modes as a function of the light 
path offset d and the aperture D 
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The turbulence profile C\ determines the correlation functions C/(6) used 
in the above equations. If C„ is unknown, higher order correlations can be 
approximately derived from the measured tip/tilt correlations in pure Kol- 
mogorov turbulence. 

Inserting C ; (6) into Eq. |5] yields the wavefront variance Cq N for N cor- 
rected modes. The isoplanatic angle 60 can then be calculated such that Oq q 1 - 
1. Fig. U shows the Strehl ratio for the wavefront variance given by Gg t + a; 
as a function of the number of corrected modes and the angle between science 
object and guide star. Because of the fast decorrelation of the higher modes, a 
low order correction can yield higher Strehl ratios at large distances between 
object and guide star (see Fig. 0). 



N ' 

2 



Strehl 
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10 corrected modes 
30 corrected modes 
100 corrected modes 



5 10 15 20 25 30 

Angle between Object and Guide Star [arcsec] 

Oq as a function of the number of 
corrected modes and the angle between science object and guide star (D = 3.5 m, single 
turbulent layer at 5000 m, ro = 60 cm) 



Figure 2. Strehl ratio for the wavefront variance o| t 



2.2. Angular Anisoplanatism of MCAO-systems 



In the case of an MCAO system with M DMs at conjugated altitudes zu the 
correlation C y (6) becomes 



/ M fc+Zi+l)/2 



cm 



£ f Cj(\zi-z\Q)C 2 n (z)dz\ / Jc 2 n {z)dz 



(9) 
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where (zo+Zi)/2 = (the integration starts from the ground) i.e. zo ■= —zu 
and Zm + Zm+i equals the upper turbulence limit of the atmosphere. For sim- 
plification it has been assumed that the number of corrected modes per area is 
the same in all layers, i.e. DMs correcting the high altitude turbulence have to 
correct more modes (because of the larger covered area) than the low altitude 
DMs. Due to the small light path offset between object and guide star in an 
MCAO system, the correlation functions have larger values than those of AO 
systems (see Fig. If there are more turbulent layers than DMs Eq. ||can be 
used to define the conjugate altitudes Zi of the individual DMs by optimizing 
Cj(Q). Inserting C/(8) into Eq. |5| yields the angular anisoplanatism o\ of an 
MCAO system (see Fig. |j). If tip/tilt is measured by using an NGS Eq. [7] has 
to be applied. 

MCAO (2 DMs) 



LGS object 




Figure 3. Anisoplanatism in an MCAO system with strong boundary layer turbulence: the 
shaded areas illustrate the amount of uncontrolled turbulence between the DMs (darker 
shading = stronger turbulence). The width of the triangle equals the light path offset 
d(z) = \zdm — z| 6- The anisoplanatism can be calculated using Eq. (9), by integrating along the 
light path the correlation functions of individual modes cj(d(z)) weighted by the turbulence 
profile C 2 n {z). 
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Figure 4. Strehl ratio for the wavefront variance o| t + Gg as a function of the angle between 
object and central LGS for a 2 DM MCAO system and LGS geometries shown in Fig. ^ (see 
section [2j|). For comparison the Strehl ratio is also given for a conventional AO system. The 
parameters are: telescope aperture D = 3.5 m, Calar Alto C^-profile, K-Band and 37 corrected 
modes. The Strehl ratio depends both on the azimuth angle a and the distance to the central 
LGS. The tip/tilt of each LGS is assumed to be known. The superiority of an MCAO system 
over an AO system, especially when seven LGSs are used, is obvious. The diameter of the 
corrected field of view can be as large as 3 arcmin. Due to the corrected cone effect, the 
MCAO system shows better performance even at the field center. 

2.3. Laser Guide Star Geometry 

To measure the wavefront distortion, several guide stars are necessary to 
cover the field. Since it is unlikely to find natural guide stars at the desired 
positions, LGSs will be necessary in most cases. Mapping the wavefront 
distortion over the desired field requires the LGSs to be pointed such that 
their light cones still overlap at the altitude of the highest turbulent layer (see 



Usually the science object does not coincide with an LGS, so different 
parts of the wavefront intersect different LGS light cones and are corrected 
accordingly. Therefore the anisoplanatism depends on the angular distance 
to different LGSs. In this paper, the fraction of intersection of each LGS 
is taken as a weighting factor for its contribution to anisoplanatism. This 
approximation leads to a linear Strehl ratio decrease close to the field center, 
see Fig. f| (instead of a more gaussian-like decrease, see Fig. |2|). However, 
it is also possible to achieve an almost flat Strehl ratio over the FOV if the 



Fig- i). 
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peak Strehl ratio at the field center is sacrificed for achieving a higher Strehl 
ratio at the edge of the field. Fig. || shows the Strehl ratio as a function of the 
azimuth angle a and the angular distance between the object and the central 
LGS. It is obvious that seven LGSs lead to a larger and more evenly covered 
field of view than five LGSs, although the Strehl ratio drop-off due to residual 
angular anisoplanatism is slightly higher. 




Figure 5. Top view at the highest turbulent layer. Five (left) and seven (right) LGSs have 
been used for the following calculations. Each circle denotes the illuminated area of the light 
cone from an LGS in the turbulent layer. If e.g. the Calar Alto turbulence profile is taken the 
Sodium layer at 90 km and the turbulent layer at 7 km result in a illuminated circle with a 
diameter of 92% of the aperture of the telescope. These parameters are used throughout this 
paper. For measuring the wavefront distortion over the desired field of view, the light cones, 
i.e. the corresponding subapertures of the wavefront sensors have to cover the field of view 
completely. The field angles used for Fig. ^ (0°, 30°, 45°) are also shown. 

The positions of the LGSs correspond to the highest accuracy for wave- 
front sensing (as in AO-systems, no angular anisoplanatism). At low altitudes, 
however, the light cones of the LGSs overlap almost completely. The wave- 
front reconstructions of the low layer are slightly different for the individual 
LGSs because the extended turbulent layers are reduced to thin layers. Since 
one has to choose one reconstruction to steer the mirror we used the central 
LGS, which defines the center of the corrected FOV. The angular anisopla- 
natism of the other LGS positions is mainly due to this low altitude correction. 
For high altitudes and little overlap between neighbouring light cones, more 
degrees of freedom allow a correction close to the optimum, as can be seen 
in Fig. ^. A more even correction can be achieved when the low altitude DM 
corrects the low altitude average wavefront error of all LGS directions. This, 
however, results in a slight loss of Strehl ratio on-axis. For this paper, we 
chose the optimal correction of the field center. 

The Strehl ratio for the wavefront variance a| t + Gg of an MCAO system 
with 2 DMs is shown as a function of the angle between science object and 
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LGS 1 LGS 2 LGS 3 




SH 3 SH 2 SH 1 



Figure 6. Side view of an MCAO geometry. Each circle in the turbulent layers indicates one 
Shack-Hartmann subaperture i.e. one actuator of the DM that is situated in the conjugate plane. 
Solid and dashed circles denote the correction for the central LGS and the other LGSs, respec- 
tively. At low altitudes, off-axis light cones mostly coincide with that of the (optimal corrected) 
central LGS, whereas at high altitudes, an independent correction of off-axis wavefronts is 
possible. 

guide star for different numbers of corrected modes (Fig. f7|) and different r$ 
(Fig. ||) . It is obvious that a high order correction at large angles 6 results in a 
higher gain for the Strehl ratio in MCAO than in single DM AO systems. For 
our calculations we assumed that the absolute tip/tilt of the LGSs is known 
(see end of section 2). Due to the (D/rn)( 5 ' 3 )-law (see Eq. ||), both AO and 
MCAO systems heavily depend on good seeing to reach a high image quality. 

Other possibilities of calculating the isoplanatic angle of MCAO systems 
are given by Tokovinin et al. (2000) and Wallner (1994). Their solutions, 
however, do not take the number of corrected modes into account. 

A number of methods have been proposed to deal with the tip/tilt de- 
termination problem. Although it is possible to use only one natural guide 
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Figure 7. Strehl ratio for the wavefront variance cjL + Og of an MCAO (2 DMs) as a function 
of the number of corrected modes and the angle between science object and optical axis 
(D = 3.5 m, 37 corrected modes, Calar Alto -profile, K-Band, 7 LGS, a = 30° (see Fig.|)). 

star (NGS) for absolute tip/tilt determination over the whole field of view 
(Johnston and Welsh, 1994), the wavefront reconstruction would not be as 
accurate as if the absolute tip/tilt was determined for each LGS. Since it is 
unlikely to find one NGS close to each LGS, the absolute tip/tilt of each LGS 
should be measured independently. Ragazzoni (1996) proposed to position 
and point auxiliary telescopes in such a way that one LGS and one NGS are 
in a line with one auxiliary telescope. Then the tip/tilt of each LGS can be 
determined by subtraction^. Another possibility is the use of polychromatic 
LGSs (Foy et al., 1995), which excite two wavelengths in the sodium layer. 
The wavelength difference of the two colours leads to a differential tip/tilt 
caused by atmospheric dispersion, from which the absolute tip/tilt can be 
calculated. Unfortunately, for a sufficient excitation of the second colour, the 
laser output power has to be increased by about two orders of magnitude. 



Since one auxiliary telescope can measure tip or tilt, two auxiliary telescopes per LGS 
are required for absolute tip/tilt-determination. 
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Figure 8. Strehl ratio for the wavefront variance o| t + Og of an MCAO (2 DMs) as a function 
of ro (scaled Calar Alto C\ -profile) and the angle between science object and optical axis 
(D = 3.5 m, 37 corrected modes, 7 LGS, a = 30° (see Fig. |)). 

3. Separating the wavefront distortion of high and low altitude 

turbulent layers 

It has been shown in the previous section that 2-DM-MCAO systems al- 
ready produce a wide and relatively evenly corrected field of view. In this 
section, we discuss methods to distinguish the wavefront aberrations of the 
corresponding two layers, a high and a low altitude layer. 

One possibility is the use of tomographic methods, as proposed by Beck- 
ers (1988). Since a large number of LGSs is required, this method is quite 
demanding for night telescopes, but can be very useful for solar telescopes 
where no LGSs are needed. New aspects of zonal tomography can be found 
in Tallon and Foy (1990). The possibility of modal tomography has recently 
been shown by Ragazzoni (Ragazzoni et al., 1999)(Ragazzoni et al., 2000). 
Solar AO systems were discussed by Scharmer et al. (2000) and Rimmele 
etal. (1999). 

In order to reduce the number of LGSsf] for measuring the individual 
wavefront distortion of two turbulent layers, the intensity information pro- 
vided in each lenslet of a SH sensor can be used, as proposed by Ribak and 
Gershnik (1996). If the SH sensor is situated in the conjugate pupil plane, 

3 However, the desired FOV should be completely covered by the LGS light cones. 
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the intensity variation at low light levels is dominated by photon noise. This 
situation can be improved by varying the position of the SH sensor as will 
be discussed in the following. For further SNR improvement and for reduc- 
ing the error caused by the wavefront modes undetectable by scintillation, a 
second SH sensor can be used, effectively forming the SH-Curvature sensor. 

Further aspects of wavefront separation and reconstruction can be found 
in (Roddier, 1988; Hickson and Burley, 1994; Ribak, 1996; Schwartz et al., 
1994). 

3.1. Introduction to scintillation 

Intensity fluctuations of star images (scintillation) are caused by the curvature 
of turbulent layers (second derivative of the phase, lensing effect, as shown in 
Fig.|). 



plane wavefront 



turbulent layer 



distorted wavefront 



aperture plane 




case 1 case 2 



focus 




Figure 9. A turbulent layer focuses (case 1) or defocuses (case 2) the incoming light, resulting 
in a larger or smaller effective aperture in the focal plane. 

Usually, astronomical observations are not affected by scintillation due 
to long integration times and large apertures. However, if intensity fluctu- 
ations of highly time-resolved measurements are not treated as noise, they 
can be used to provide information about the altitude distribution of turbu- 
lent layers. Good introductions to the theory of scintillation can be found in 
(Dravins et al., 1997a; Dravins et al., 1997b; Dravins et al., 1998; Elsasser 
and Siedentopf, 1959; Reiger, 1963). 
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The variance of the measured intensity a 2 consists of the variance of the 
scintillation a| and other noise sources (detector noise and photon noise 



aj = a 2 s + a 2 D + o 2 P . (10) 



If the detector noise and the photon noise are known, the variance of the 
scintillation can be determined. According to Reiger (1963), the theoretical 
value of a| is given by 

ajocD- 7 / 3 (cosy)- 3 J z 2 C 2 n (z)dz, (11) 

D and y denoting the telescope aperture and the zenith angle, respectively. 
Simulations with the program Turbulence (Berkefeld et al., ) have shown that 
the normalized variance of the intensity for a single layer at altitude z can be 
approximated by 

a 2 = D-y\ 2 z 2 r 5/3 (cosy)- 3 . (12) 

Because of ro oc X 6 / 5 , the scintillation does not depend on the wavelength. For 
an extended turbulence profile one obtains 

oj = 16.7-D- 7 / 3 (cosy)- 3 y z 2 C 2 n (z)dz, (13) 

Due to the factor z 2 the scintillation is mostly caused by high altitude turbulent 
layers. For 0$ > 10%, e.g. for large zenith angles y or for small telescope aper- 
tures D, the scintillation begins to become nonlinear, approaching a maximum 



value, and Eq. [13| is no longer valid (Dravins et al., 1997a; Protheroe, 1954). 



3.2. Wavefront sensors for separating the contributions of 
high and low altitude turbulent layers 

Since the wavefront gradient in each SH subaperture is the sum of the wave- 
front distortions of all turbulent layers, additional information is needed for 
distinguishing the influence of individual layers. By using the spatially re- 
solved scintillation provided in each subaperture of the SH sensor, the distor- 
tion of one of the layers can be reconstructed. Together with the known sum 
of the aberrations, the wavefront distortion of the other layer can then also be 



determined (Fig. 10). 



As in a conventional AO-setup, the SH sensor measures the tip/tilt in each 
subaperture, thus delivering the integrated wavefront error. Because the SH 
sensor is situated in the conjugate plane of one of the turbulent layers (sub- 
sequently named first layer), this layer has no effect on the intensity of the 
SH pattern. The - for the SH sensor - defocused other (subsequently named 
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Shack-Hartmann sensors 
in the conjugate plane of one 
of the turbulent layers 

Figure 10. Setup of a wavefront sensor separating the wavefront distortions of two turbulent 
layers. For clarity, the light paths of only two guide stars have been plotted. One SH sensor 
per LGS is required. 



second) layer leads to intensity fluctuations /^(x^) i n tne subapertures from 
which the wavefront distortion (j) of the second layer can be reconstructed 
(Ribak and Gershnik, 1996; Glindemann and Berkefeld, 1996): 

where indices denote the layers, z is the altitude coordinate, k = Itz/X, and 
dashed variables indicate the conjugate planes. The radial wavefront tilt af- 
fects the intensity measurement at the edge of the aperture (Roddier, 1988). 
Although the edge radial tilts can be extracted directly from the SH tilt mea- 
surement and thus their contribution can easily be corrected, they must be 
known for each layer. They can be estimated by first determining the mean 
gradient m; (averaged over all LGSs j) of each edge subaperture i. The m, 
can be regarded as the edge gradients of the low altitude turbulence (which 
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are the same for all LGSs), whereas its deviation r,j can be regarded as the 
high altitude gradients 

n 

r U = l U - m < with m i = ■ > ( 15 ) 

n 

with ttj being the measured gradients. It is obvious that this is only an ap- 
proximation which will be more accurate with a higher number of LGS (and 
thus a larger corrected field of view). Once the wavefront distortion of one 
of the layers has been reconstructed, subtraction from the measured sum of 
the wavefront errors delivers the wavefront distortion of the other layer. The 
error made by averaging the radial gradients, as described by Eq. is still 
unknown. This problem should therefore be addressed more closely in the 
future. 

One can define the contribution of an individual layer i to the scintillation 
in a conjugate plane of altitude h by (see Eq. |T|) 

a 2 Sii (h) = 16.7.D s - u 7 b /3 (cosYr 3 / {z-h) 2 C 2 t {z)dz, (16) 

Li 

with D su b the subaperture size of the SH sensor that is placed in the conjugate 
plane of altitude h. 

The goal of the measurement is to determine the contribution of the more 
turbulent of the two layers. Thus, the SH sensor is placed in a conjugate plane 
where the scintillation of the more turbulent layer (in this case layer 2) is very 
large and the contribution of layer 1 is very small. Then, one can write the 
signal-to-noise ratio (SNR) as the normalized mean square error (1/SNR 2 ) of 
the intensity measurement as: 

i 

1=— • (17) 

SNR 2 a 2 s2 (h) 

The numerator consists of the photon noise \/N 1 / 2 , the read noise R ■ n l l 2 
(n is the number of pixels used for the intensity measurements, R denotes the 
readout noise per pixel), and the contribution of layer 1 that is treated as noise 
for the measurement of layer 2 that constitutes the signal in the denominator. 
The scintillation error due to the non discrete layering of the C 2 -profile is very 
small and need not be taken into account (Berkefeld, 1998). 

Typically the measured scintillation is of the order of a few percent. This 
poses problems at low light levels where the scintillation competes with the 
shot noise. Although it is possible to measure the scintillation of the more 
turbulent layer (usually the ground layer) and thus get a smaller relative error, 
this does not lead to a better separation of the layers because the absolute 
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error of the wavefront reconstruction is also proportional to the turbulence 
strength. Therefore it does not matter whether the scintillation of the strong 
or the weak turbulence is measured. Instead, the mean square error can be 
reduced by moving the SH sensor further away from the conjugate plane of 
the first turbulent layer, as shown in Fig. 1 1 . 

MSE 

1 r 

/ 

/ 

/ 

0.8- „ ^ ^ 




1 = 7 mag 

-5 -4 -3 -2-10 1 

conjugate height of Shack-Hartmann sensor [km] 

Figure 11. Normalized mean square error (1/SNR 2 ) of the scintillation measurements as a 
function of the conjugate altitude of the SH sensor for different LGS intensities. The sub- 
aperture size and the read noise have been assumed as D su b = 0.5 m (35-40 corrected modes) 
n = 16 pixels and R = 3e~, which are typical values for the AO system ALFA at the Calar 
Alto 3.5 m telescope. The integration time of the wavefront sensor was chosen to be optimal 
for each LGS brightness. In this figure, the scintillation caused by the high altitude turbulence 
is measured. Increasing LGS brightness leads to increasing optimal altitude of the SH sensor. 
In order to have an SNR better than 4 the mean square error has to be smaller than 1/16. 



Although this leads to an unwanted contribution to scintillation from the 
now defocused first layer, the SNR increases within certain limits due to the 
much higher scintillation signal from the second layer. However, in the case 
of a 589 nm LGS with a 10 mag G-star brightness equivalent, the maximum 
SNR is 1.15, which is not sufficient for wavefront separation. 

3.3. The Shack-Hartmann-Curvature sensor 

For low light levels the amount of scintillation has to be increased by mov- 
ing the SH sensor further away from the conjugate planes of the turbulent 
layers. Since the scintillation effects of two layers cannot be distinguished 
with one SH sensor, an additional SH sensor is necessary, forming the Shack- 
Hartmann-Curvature sensor. In order to refine the SHC setup presented by 
Glindemann and Berkefeld (1996) we propose to position the individual SH 
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sensors in such a way that a mean intensity fluctuation of 10% (the non- 
linearity limit of scintillation) is achieved in each SH sensor. This is the case 
at the conjugate planes of a very large positive altitude and a large negative 
altitude, leading to an excess of illumination in one plane and to a lack of 
illumination in the other, similar to the Curvature Sensor (Roddier, 1988). 
For the Calar Alto C\ -profile this results in altitudes of approx. 16 km and 
-12 km. 

MSE 

1 



0.8 - 



0.6- 



. 4 



0.2 



7 8 9 10 11 12 

Guide Star Brightness [mag] 

Figure 12. Normalized mean square error of the intensity measurements performed by a SHC 
sensor as a function of the guide star brightness. The conjugate altitudes of the SH sensors 
have been assumed as 1 6 km and - 1 2 km, the other parameters are the same as in fig 

The scintillation variance a| shc(^i > ^2) of a SHC sensor can be expressed 
as the sum of two SH sensor signals each of which receives contributions of 
high and low altitude turbulence. For turbulence separation, the signal of each 
SH sensor is the difference of high and low altitude turbulence scintillation: 



0|,SHC(*1,*2) = 16 - 7 « 

(j\( z - hl f + ( z -h 2 f)C 2 n (z)dz - 

L, L 2 



-7/3 
sub 



(cosy) 



-3 



(18) 



{{z-h l f + (z-h 2 f)cl{z)dz) 



Then, the mean square error for a SHC is 
1/SNR 2 = 



_L 4_ n R 



(19) 



°S,SHC(^1^2) 

An LGS brightness of 8 mag leads to an SNR of 4 which should be sufficient 



for most applications (see Fig. 12). Therefore observations in the NIR and 
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a typical LGS brightness of 10 mag require one SHC sensor per LGS. Ob- 
servations in the visible require much more powerful LGSs in order to sense 
and correct a higher number of modes. Then, a single SH per LGS can be 
used. In any case, the -profile and the zenith angle should be monitored 
continuously in order to adjust the SH sensors and the DMs accordingly. 

There is another aspect that makes the use of a SHC sensor more desirable 
than a SH sensor: Since the scintillation signals of the wavefront modes are 
not linearly independent, it is not possible to measure the wavefront modes of 
zero Laplacian with a SH sensor directly. Although it is possible to estimate 
these modes by averaging the radial wavefront tilts as described in section 
3.2 , the remaining error should be reduced further by using the intensity 
measurements for applying a priori (statistical) knowledge. Since a SHC sen- 
sor provides two and a SH sensor only one measurements per subaperture, 
the application of a priori knowledge will work better with a SHC sensor. 
However, for deriving a stable algorithm that operates in real-time, more de- 
tailed work has to be done in the field of wavefront reconstruction by intensity 
measurements. 



4. Conclusion 

By using MCAO, it is possible to overcome the most severe disadvantage 
of AO, the very small corrected FOV. We have shown how the geometry of 
MCAO systems affects the angular anisoplanatism and thus the size of the 
corrected FOV. A setup with seven LGSs for wavefront sensing leads to a 
wide and relatively evenly corrected field. In the case of the Calar Alto 3.5 m 
telescope, this setup would lead to a FOV of about three arcminutes. Eq. || 
plays an important role in calculating the wavefront error caused by angular 
anisoplanatism. By its minimization one obtains the optimal position of the 
DMs (according to the C% -profile). Since the remaining anisoplanatism inside 
the field is rather small, a wavefront correction with two DMs seems to be a 
good compromise between anisoplanatism, cost and system complexity. The 
separation of the wavefront errors of the two layers can be accomplished in 
various ways: For solar telescopes, tomographic methods, as proposed by 
Beckers, seem to be the most accurate and easy to implement way. For night 
telescopes, the separation can be done by using the intensity information 
provided by each lenslet of the SH sensor or SHC sensor. The SHC sensor 
should be preferred because it allows a fainter limiting magnitude and can 
reduce the separation error caused by the wavefront modes with zero Laplace 
operator. 

The high costs for the LGS setup and the absolute tip/tilt determination 
will prevent MCAO systems from being used at existing 3.5 m-class-telescopes, 
at least as long as other improvements of conventional AO systems are possi- 
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ble. At 8+ m-class-telescopes however, focal anisoplanatism will require the 
use of multiple laser guide stars. Furthermore, the cost of MCAO systems 
compared to those of the telescope will decrease, so that MCAO will become 
a common feature at large telescopes. 



Appendix 
A. Calar Alto C^-profile 



Fig. 13 shows a slightly simplified C\ -profile measured by Kliickers et al. 
(1998) at the Calar Alto Observatory, Spain. The upper turbulent layer at an 
altitude of 7 km delivers the main contribution to the angular anisoplanatism, 
the lower turbulent layer determines most of r$. For a 2-DM-MCAO, the 
optimal altitudes for the DMs are 400 m and 6900 m. It should be noted, that 
the turbulence profile can change rather quickly. Kliickers reported a change 
of the upper turbulent layer strength by a factor of two in only a few minutes. 
Therefore, frequent C\ -measurements should be made to adjust the optimum 
altitude of the DMs. 



^ 2 r -2/3, 

C„[m 1 



4000 6000 

height [m] 



Figure 13. C„ -profile at the Calar Alto Observatory (Klueckers et al., 1998), ro = 66 cm in 
K-band) 
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